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The contribution of the light-by-light diagram to the g factor of electron and muon bound in 
Coulomb field is obtained. For electron in a ground state, our results are in good agreement with 
the results of other authors obtained numerically for large Z. For relatively small Z our results have 
essentially higher accuracy as compared to the previous ones. For muonic atoms, the contribution 
' , is obtained for the first time with the high accuracy in whole region of Z. 
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The progress in experimental investigations of the g factor of a bound electron (Ij and muon y| in ions stimulated 

' intensive theoretical investigation of various contributions to this quantity. The contributions of self-energy, vacuum 

K>" polarization, and nuclear effects have been considered 0, 0, S S S An essential part of the theoretical 

' uncertainty has been related to the contribution of the vacuum polarization of an external homogeneous magnetic 

i—j I field in the electric field of atom (so-called the "magnetic-loop" contribution) . The corresponding diagram is shown 

. in Fig. n In this diagram, double line in the fermion loop corresponds to the electron propagator in the Coulomb 

I field. Note that the contribution of the free electron loop to the vacuum polarization of a homogeneous magnetic 

. field vanishes due to the gauge invariance. The first non-vanishing terms of expansion with respect to the Coulomb 

CD ' field shown in Fig. ^ is the contribution of virtual light-by-light scattering with one of the quanta corresponding 

, to the external magnetic field. The results of numerical calculations of the magnetic-loop contribution, which take 

O • into account all orders of the parameter Za (Z is the nuclear charge number, a = is the fine-structure constant, 

I \ S = c = 1), are presented in Ref. 0. At present, the most accurate experimental data are obtained in the region 
of medium Z. Unfortunately, in this region the uncertainty of the results of Ref. is very big, being, e.g., 100% 

(— I for Z = 12. In Ref. 9], the leading in Za magnetic- loop contribution to the g factor of an electron in S state of a 

• • . hydrogen-like ion has been derived. It reads 
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I. INTRODUCTION 
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where go is the Lande factor equal to two for S state. One can compare this correction with the result of [6j for 
rather large Z where the accuracy of the numerical calculation is reasonable. This comparison shows the noticeable 
difference which can be attributed to the contribution of the next-to-leading terms in Za-expansion, starting from 
aiZaf . Since the numerical factor in Eq. |^ is very small (~ 1/30), the next-to-leading terms could give a noticeable 
contribution to the g factor even at small Z, if the corresponding numerical factor is of order of unity. 

In the present paper, we generalize Eq. to the case of arbitrary bound electron state. We also calculate the 
next-to-leading contribution of magnetic loop to the g factor of the electron in arbitrary state (or the magnetic moment 
of the electron in this state). It has the form Agi ~ a(ZoLf{a\ ln(l/ZQ;) + a-i), where a\^2 are some constants and a\ 
is not zero only for S states. In order to calculate this contribution, it is sufficient to take into account the diagrams 
of virtual light-by-light scattering and use the nonrelativistic wave functions of the bound electron. Comparison of 
the correction Ago + Agi for XS^ji state with the results of '3| shows that the account of Agi does not provide good 
agreement for relatively small Z ~ 30, where the numerical calculations were performed with sufficient accuracy. 
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Thus, for such Z it is necessary to take into account next terms in Za.. These terms have two different origins. 
First, they come from the relativistic corrections to the wave function of a bound electron. Next, they come from the 
higher-order contributions to the electron loop. Note that the diagram in Fig. QJcan be interpreted as the contribution 
of the scattering of the magnetic quantum in a Coulomb field (virtual Delbriick scattering) to the g factor. It is known 
that the Coulomb corrections to the Delbriick amplitude for momentum of quantum q "^m (m is the electron mass) 
are numerically small even for large Z [TlIIT2i|. In contrast, the account of the corrections to the wave function is very 
important, starting from relatively small Z . We calculate the correction /S.g using the relativistic wave function and 
the leading approximation for the electron loop. As a result we have obtained good agreement with the numerical data 
of 1^ even for very large Z (difference is 4% for Z = 92). Using such approach, we have calculated the corresponding 
correction of the electron loop to the g factor of a bound muon. 




- X + • • • 



FIG. 1: The diagram corresponding to the magnetic- loop contribution to the g factor of a bound electron and first non- vanishing 
terms of expansion of this loop with respect to the Coulomb field. Double line denotes the electron propagator and the wave 
function in a Coulomb field, the dashed line with the cross denotes the Coulomb field, the wavy line with the square denotes 
the external homogeneous magnetic field, the internal wavy line corresponds to the photon propagator. 



II. GENERAL RELATIONS 



Let us consider the amplitude T of interaction of homogeneous magnetic field B with the electron bound in a 
hydrogenlike ion. In the zero approximation, it reads (see, e.g., |l3l |) 



T''' - e I ■ jZ = ^^i^ / drr^Mr)Mr) , (2) 



dk _ ekB ■ (J ) 

(2^ J(J + 1 

where jk is the Fourier transform of the electron current j{r) = ilj{r)^^{r), the wave function ij; has the form 

n is the spherical spinor |14| with the angular momentum J and orbital momentum L, fl — —(cr ■ n)fl, k — {J + 
l/2)sign(i — J). In Eq. jSJ we have used the relation 

ik X Ak = {2TT)^Sik)B . (4) 

Note that a sign of T^^^ is opposite to that of Hamiltonian. Substituting the radial wave functions fi{r) and f2ir) 
for the Coulomb field (see, e.g., ^^), we obtain for the arbitrary bound state 

2 m 

2k, ( 2Ke\ 2k f 2k \ 
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\ m J f - 4kM ^l + iZaY/i-f + nrf 



where rir is the radial quantum number, e is the binding energy, and 7 — ^ k? — (Za)'^. The particular cases of this 
formula obtained earlier are presented in Ts^. In the non-relativistic approximation {Za ^ 1), Eq. (0) turns to 

(0) eB ■ (J) 2k 

We now pass to the calculation of the amplitude T^^^ corresponding to the diagram shown in Fig. ^ It has the 
form 

J (2^)3 J (2^)3 q2 ^k^^l Jq , [') 



3 



where the amphtude A4 oi the virtual Delbriick scattering in the case k <^m has the form following from the gauge 
invariance 

i (J, • q) - q'k^] T{q/m, Za) . (8) 

TO-' 

Note that T is even function of Za. In the leading in Za approximation (contribution of light-by-light scattering), 

T{q/m, Za) = {ZafF{q/ra) , (9) 
with F{Q) = 7/1152, see Ref. 0. From Eqs. Q, ©, ©, and © we obtain 

OO OO 

= /^^(^^)-^ - (J) h,F[q/m) [drrMr)Mr) ( - cos,.) (10) 
■Km-^J{J +1) J J \ qr ) 



Using the relation A^" = [2a{Za)'^{k ■ q)/m^]F{q/m), following from Eq. ©, and the gauge invariance of the 
light-by-light scattering amplitude, we can represent F{q/m) in the following form 



F{qM = -j 



dQ q-(VfcM)|fe=o 
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27r J Q^{q-QY q 
M = 2tJ ^Sp|g(p)7'G(p-A:)7° 



G(p + Q-g)f G(p-t-g)7° 



+G{p + Q- g)7°G(p - g)y + G(p - Q - fc)7"G(p - 9)7' 



(11) 



where G{p) = [p — m]~^ is a free electron propagator. Straightforward calculation leads to the representation of the 
function F in the form of two-fold integral with respect to the Feynman parameters. Resulting formulas being rather 
cumbersome are not presented here explicitly. For x — q/m <^ 1, the first two terms of expansion of the function 
F{x) have the form 

F(x) = ^(1 + —x). (12) 
^ ' 1152^ 35 ^ ' 

The first term in this formula agrees with the result of Ref. 9]. For x ^ 1, the asymptotics of the function F{x) 
reads 

Fix) = ^ . (13) 

For arbitrary x, we performed the numerical tabulation of the function F{x). The result is shown in Fig. |21and in 
Table D 



III. CORRECTION TO g FACTOR AT SMALL Za 

In order to obtain the leading term of expansion in Za of the amplitude T^-^\ it is sufficient to use Eq. I|1U|) 
with the substitution F(q/m) — > F(0) and the wave functions taken in the nonrelativistic approximation. In this 
approximation /i(r) coincides with R{r), the radial part of the nonrelativistic wave function, and 

f2ir)^^(R'ir) + ^R{r)) . (14) 



(15) 



2m \ r 

The correction Ag to the Landc factor is determined by the relation 

Ag _ r(i) 

50 " r„(°) ■ 



Taking in Ea. ()10|l the integral over g, and then over r, we obtain the leading contribution Ago for the arbitrary 
state 

Ago ^ 7a(Za)5 ^ lajZa)^ 

go 144n3(2L-Hl)K(2K- 1) 288^3 J( J + i)(2J + 1) ' ^ ' 
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FIG. 2: The ratio F{x)/F{0) as a function of a; = g/m 



where n — rir + \k\ is a principal quantum number. For S states {L ~ 0, k = —1), this resuh is in agreement with 
Eq. (|TJ obtained in Ref. ^j. 

The relativistic corrections to the wave function as well as the corrections to the magnetic loop have the relative 
magnitude {Za)^. Therefore, the term Agi of the order a(Za)^ can also be obtained with the use of the nonrelativistic 
wave functions and magnetic loop in the leading approximation (light-by- light scattering diagrams). For L 7^ 0, it is 
sufficient to substitute the second term of expansion of F{x), see Eq. (|12l) . in Eq. (|10l) . Then we obtain 



Agi _ 2a{Zaf 

~ 457rn3(2L-f 1)(2k- 1)2 \L{L + 1) 



1 



(17) 



For S states, calculation of A51 is more complicated. For nS state, /i(r)/2(r) = {n/m)p'^{r), where p„(r) is the 
electron density in the nonrelativistic approximation. Substitution of Eq. H12|) in Eq. (|10|l leads to logarithmic diver- 
gence. Therefore, it is convenient to split the region of integration over r in Eq. p(J|l into two: [0, tq] and [ro, 00) with 
1/m <^ ro l/{mZa). In the first region, we can replace p'{r) by p'(0) and take the integral over r. In the second 
region, we can use the expansion H12I) and take the integral over q. Sum of these two contributions, as it should be, 
is independent of tq. The final result reads 



Agi 
90 



MZaff 1 
r "^7; a -On I , 



1357rn'' 



Za 



00 

-- + ^ / dx\nxF"{x) ^2.6, 
2 8 J 



bn 



-C 



P'rM 



dr \n{mZar)pn{r) . 



(18) 



where C — 0.577... is the Euler constant. For each n, the coefficient &„ can be easily calculated so that bi = ln2 
0.693, &2 = 5/8 = 0.625, 63 55/54 + In 2/3 « 0.613, 600 = C + In 2 - 2/3 « 0.604. 



IV. CORRECTION TO g FACTOR AT Za ~ 1 

As it was pointed out in the Introduction, the sum A^q + Agi gives a good approximation to Ag only for small 
Z. For intermediate Z, it is necessary to account for the next terms in Za. The largest corrections are due to the 
significant difference between the relativistic wave function and the nonrelativistic one already at intermediate Z. At 
the same time, the difference between the function and its leading approximation {Za)'^F results in the corrections 
which are numerically small even for large Z. Using the numerical results for F(x) and the relativistic wave functions, 
we have performed the tabulation of Ag for various Z, using T^^^ from Eq. H10|l as an approximation to T*-^). The 
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results of this tabulation for IS1/2, 251/2, and 2Pi/2 states are presented in Table UTI For lS'1/2 state, we also present 
the contribution of the first two terms of expansion in Za, Eqs. H18I) . (|16|1 . and the correction Ag^r obtained with 
the use of nonrelativistic wave functions. The results for lS'1/2 are also shown in Fig|3| 

.04, ■ 




. 
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FIG. 3: The correction Ag for lS'1/2. Solid curve: the exact result; dashed curve: Ago + Agi; dotted curve: A(;„,.. 

For Z < 10, both Ago + Agi and Agnr coincide with Ag with accuracy better than one percent. The difference 
grows with Z, reaching 10% at Z '--^ 30 for Ago + Agi and at Z ^ 50 for Apnr- 

In Table Hll we also show the results of numerical tabulation from Ref. p\ for 151/2 state. For 30 < Z < 70, our 
result for agrees with that obtained in Ref. ^] within 1 2% percent. The difference between these two results for 
Z < 30 is due to poor accuracy of the numerical results of Ref. [a| . For Z > 70 the difference increases and becomes 
8% for Z = 92. This difference corresponds to the contribution of next-to-leading terms in magnetic loop, that was 
taken into account in Ref. 6] and omitted in our paper. Thus, the effect of these terms is small in a wide region of 
Z, while the relativistic effects in the wave function become important already at relatively small Z. 



V. THE CORRECTION Ag FOR MUONIC ATOMS 

The correction Ag to the g factor of a bound muon due to the electron magnetic loop can be obtained from Eq. H10() 
with fi{r) and /2(f) being the wave functions of the muon. The asymptotics of Ag for ^Za/ (mn^) « 1.5Z/n? ^ 1 (/i 
is the muon mass) can be calculated as follows. We split the region of integration over q in Eq. 1)1 0|l into two: [0, qo] 
and [qo, 00) with m <C go ^ iiZa/n"^. In the first region we can replace [{qr)~^ singr — cosqr] by ((jr)^/3 and take 
the integral over r. In the second region we can use the asymptotics Eq. (|13|l and take the integral over q. Summing 
these two contributions, we obtain 

Affa. = g ^"^5"^' [ln(/iZa/m) - A - B] , 

A = 2 dy\nydy{y''F{y)) ^2.24, 
Jo 

4 4 ^ ^ 

^ = ""-3- Zail-2ns/,) i Inx , 

/i(x) = i^lZa)-''^^f^{x/^lZa), f2ix)^i^lZa)-^/^f2ix/^,Za), (19) 
g is defined in Eq. (0. For 15i/2 state we obtain 

5 = ^(1 + 27), B = C-^+V(27 + 2)-ln2, j ^ y/l^{Zaf . (20) 
For n = Ur + \k\ 1, we have 

2k 

ff-.9o = :^— 7, B = C + Hn'/2). (21) 

The formula H19|l can be interpreted as follows. In Ref. the logarithmic contribution of the electron vacuum 
polarization to the magnetic moment of a heavy nucleus was calculated. The result obtained has the form 

Ag 2a(Za)'^H(Za) , , , , 

— = ^ 5 ^ ^ ln(l/mi?„„,i) , (22) 

g Stt 
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where Rnuci is the nuclear radius, Rnuci <IC l/m. The coefScient {Za)'^H(Za) was calculated exactly in Za, i.e., with 
the account of all Coulomb corrections to the electron loop. The function H{Za) tends to unity when Za 0, and 
significantly differs from unity only for very large Z. Large logarithm ln{l / mRnuci) in Eq. (|22|l appears as a result of 
integration over distance r in the region Rnuci J' <C 1/m. We can consider the muonic atom as some nucleus with 
the effective radius Rnuci ^ r? j [iZa. In the case ^Za/{mn'^) ^ 1, we have Rnuci ^ l/m. Substituting this radius 
into Eq. (|22|l and replacing H{Za) 1 (that corresponds to the contribution of light-by-light scattering), we obtain 
the logarithmically amplified term in Eq. H19|) . Note that the coefficient in Rnuci corresponds to the asymptotics 
of B in Eq. (|19|l at n ^ 1. Strictly speaking, the charge of such effective nucleus is Z — 1, but not Z. However, under 
the condition ^,Za/ {rmi?) « 1.5Z/n^ ^ 1, this difference is not important. 

In Table imi we present Ag for IS1/2 state of muonic atom calculated for arbitrary Z. For comparison, we present 
also asymptotics Eq. (fT^ . As it should be, the accuracy of asymptotics (fT^ increases with Z being 4% for Z = 40 
and 1% for Z = 92. 

In summary, we found higher-order magnetic loop corrections to the bound g factor in order a{Za)^ for arbitrary 
state. Despite a small coefficient in the leading term of order a{Za)^, and the logarithmic enhancement of the 
higher-order contribution, the leading term still dominate for Z — 6 and Z = 8, important for experiment. Previously 
used numerical results show a certain underestimation of the magnetic-loop contribution for Z < 20. Theoretical 
description of this contribution presented in this paper is more reliable. The difference of less than a few percent 
between our analytic results and the numerical calculations of Ref. 6] at high Z (80 ~ 90) shows that the contribution 
of the higher-order terms in magnetic loop may be safely neglected for Z < 50. We also calculated the correction 
Ag for the bound muon, and its behavior is very peculiar. All known contributions to the bound g factor scale as 
or n~^. The correction found in this paper does not contain such a strong suppression factor. This correction is 
a dominant bound state QED correction for a bound muon, which even for 15*1/2 state supersedes the free vacuum 
polarization term . The results obtained significantly diminish the uncertainty of the theoretical predictions for the 
g factor of a bound particle. 

This work was supported in part by Russian Science Support Foundation, RFBR Grants Nos. 03-02-16510 and 
03-02-16843. 
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TABLE I: Function F{x)/F{0) versus x = q/m; F{0) = 7/1152. 



9 





l'S'i/2 


251/2 


2P1/2 


z 


Ago + 


1 




A<7 


Ag(Ref. [6]) 


8A3 


8A5 


1 


4.935 X 10 


-9 


4.934 X 10 


-9 


4.934 X 10 


-9 




4.936 X 10 


-9 


1.638 X 10 


-9 


2 


1.58 X IQ- 


7 


1.58 X 10" 


-7 


1.58 X 10" 


-7 




1.58 X 10" 


-7 


5.26 X 10" 


8 


3 


1.2 X 10" 


6 


1.2 X 10" 


6 


1.2 X 10" 


6 




1.2 X 10" 


6 


4.01 X 10" 


7 


4 


5.04 X 10" 


6 


5.04 X 10' 


-6 


5.04 X 10" 


-6 




5.05 X 10" 


-6 


1.69 X 10" 


6 


5 


1.53 X IQ- 


5 


1.53 X 10' 


-5 


1.54 X 10" 


-5 




1.54 X 10" 


-5 


5.18 X 10" 


6 


6 


3.79 X IQ- 


5 


3.8 X 10" 


5 


3.81 X 10" 


-5 




3.82 X 10" 


-5 


1.29 X 10" 


5 


7 


8.16 X 10" 


5 


8.17 X 10" 


-5 


8.2 X 10" 


5 




8.23 X 10" 


-5 


2.81 X 10" 


5 


8 


1.58 X 10" 


4 


1.58 X 10" 


-4 


1.59 X 10" 


-4 




1.6 X 10" 


4 


5.49 X 10" 


5 


9 


2.83 X 10" 


4 


2.84 X 10" 


-4 


2.86 X 10" 


-4 




2.87 X 10" 


-4 


9.94 X 10" 


5 


10 


4.76 X 10" 


4 


4.78 X 10" 


-4 


4.82 X 10" 


-4 




4.84 X 10" 


-4 


1.69 X 10" 


4 


11 


7.61 X 10" 


4 


7.66 X 10" 


-4 


7.72 X 10" 


-4 


3(3) X 10"'' 


7.76 X 10" 


-4 


2.73 X 10" 


4 


12 


1.17 X 10" 


3 


1.18 X 10" 


-3 


1.19 X 10" 


-3 


4(5) X 10"'' 


1.19 X 10" 


-3 


4.24 X 10" 


4 


13 


1.72 X 10" 


3 


1.74 X 10" 


-3 


1.76 X 10" 


-3 


8(5) X 10"'' 


1.77 X 10" 


-3 


6.35 X 10" 


4 


14 


2.48 X 10" 


3 


2.51 X 10" 


-3 


2.54 X 10" 


-3 


1.4(1.0) X 10"^ 


2.56 X 10" 


-3 


9.25 X 10" 


4 


15 


3.46 X 10" 


3 


3.52 X 10" 


-3 


3.57 X 10" 


-3 


2(1) X 10"^ 


3.6 X 10" 


3 


1.31 X 10" 


3 


16 


4.74 X 10" 


3 


4.82 X 10" 


-3 


4.9 X 10" 


3 


3(1) X 10"^ 


4.95 X 10" 


-3 


1.82 X 10" 


3 


17 


6.35 X 10" 


3 


6.48 X 10" 


-3 


6.6 X 10" 


3 


5(2) X 10"^ 


6.67 X 10" 


-3 


2.48 X 10" 


3 


18 


8.36 X 10" 


3 


8.56 X 10" 


-3 


8.73 X 10" 


-3 


6(2) X 10"^ 


8.83 X 10" 


-3 


3.31 X 10" 


3 


20 


1.39 X 10" 


2 


1.43 X 10" 


-2 


1.46 X 10" 


-2 


1.0(3) X 10"^ 


1.48 X 10" 


-2 


5.66 X 10" 


3 


24 


3.28 X 10" 


2 


3.45 X 10" 


-2 


3.56 X 10" 


-2 


3.3(3) X 10-2 


3.63 X 10" 


-2 


1.44 X 10" 


2 


28 


6.72 X 10" 


2 


7.22 X 10" 


-2 


7.53 X 10" 


-2 


6.9(3) X 10"^ 


7.7 X 10" 


2 


3.19 X 10" 


2 


32 


0.123 




0.136 




0.144 




0.138 


0.148 




6.37 X 10" 


2 


36 


0.207 




0.238 




0.254 




0.249 


0.262 




0.118 




40 


0.325 




0.389 




0.421 




0.410 


0.437 




0.206 




44 


0.481 




0.607 




0.665 




0.658 


0.695 




0.341 




48 


0.676 




0.907 




1.01 




1.01 


1.06 




0.545 




52 


0.904 




1.31 




1.48 




1.48 


1.56 




0.841 




56 


1.15 




1.84 




2.1 




2.12 


2.24 




1.26 




60 


1.41 




2.51 




2.92 




2.95 


3.13 




1.85 




64 


1.63 




3.35 




3.97 




4.03 


4.29 




2.66 




68 


1.77 




4.4 




5.3 




5.39 


5.77 




3.76 




72 


1.78 




5.67 




6.96 




7.11 


7.62 




5.23 




76 


1.55 




7.2 




9. 




9.24 


9.93 




7.18 




80 


0.983 




9.02 




11.5 




11.9 


12.8 




9.75 




83 


0.252 




10.6 




13.7 




14.2 


15.3 




12.2 




88 


-1.77 




13.7 




18.1 




18.9 


20.5 




17.5 




92 


-4.34 




16.5 




22.5 




23.5 


25.5 




23.1 





TABLE 11: The quantity r? in units 10"®, calculated in various approximations for 151/2, 25i/2, and 2Pi/2 states. Our results 
are obtained with the account for the magnetic loop in the leading approximation (contribution of light-by- light scattering). 
The quantity /S.gnr denotes the correction obtained with the use of Eq. llOU with the functions /i(r) and /2(r) taken in the 
nonrelativistic approximation, see Eq. 1141 . 



z 


Aff 


Agas 


Z 


Afl 


Agas 


1 


1.043 X 10"^ 


-0.2701 


24 


159.5 


145.9 


2 


0.1274 


-0.6226 


28 


233.7 


217.9 


3 


0.484 


-0.7983 


32 


324. 


306.1 


4 


1.186 


-0.6592 


36 


430.8 


411. 


5 


2.317 


-0.109 


40 


554.4 


532.8 


6 


3.944 


0.9268 


44 


694.9 


671.5 


7 


6.124 


2.508 


48 


852.3 


827.4 


8 


8.904 


4.687 


52 


1026. 


1000. 


9 


12.33 


7.506 


56 


1217. 


1189. 


10 


16.43 


11. 


60 


1424. 


1395. 


11 


21.24 


15.22 


64 


1646. 


1616. 


12 


26.8 


20.17 


68 


1883. 


1852. 


13 


33.13 


25.9 


72 


2134. 


2103. 


14 


40.26 


32.42 


76 


2398. 


2366. 


15 


48.2 


39.77 


80 


2673. 


2641. 


16 


56.98 


47.96 


83 


2886. 


2854. 


17 


66.62 


57.01 


88 


3251. 


3219. 


18 


77.14 


66.94 


90 


3400. 


3368. 


20 


100.9 


89.52 


92 


3550. 


3519. 



TABLE III: Ag in units 10"^ for IS1/2 state of muonic atom. Agas is the asymptotics l|19|l . 



